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Abstract 

Background: Prior animal model and human-based studies have linked selenium concentrations to decreased risk 
for depression; however, this work has not focused on household groundwater levels or specific depressive 
symptoms. The current study evaluated the link between groundwater selenium levels and depression. We also 
sought to determine if a functional polymorphism in the glutathione peroxidase 1 (GPX1) gene impacted this link. 

Methods: We used a cross-sectional design to analyze data from 585 participants (183 men and 402 women) from 
Project FRONTIER, a study of rural health in West Texas. Residential selenium concentrations were estimated using 
Geospatial Information System (GIS) analyses. Linear regression models were created using Geriatric Depression 
Scale (GDS-30) total and subfactor scores as outcome variables and selenium concentrations as predictor variables. 
Analyses were re-run after stratification of the sample on GPX1 Pro198Leu genotype (rsl 050454). 

Results: Selenium levels were significantly and negatively related to all GDS and subfactor scores accounting for up 
to 17% of the variance beyond covariates. Selenium was most strongly protective against depression among 
homozygous carriers of the C allele at the Pro198Leu polymorphism of the GPX1 gene. Analyses also point towards 
a gene-environmental interaction between selenium exposure and GPX1 polymorphism. 

Conclusion: Our results support the link between groundwater selenium levels and decreased depression 
symptoms. These findings also highlight the need to consider the genetics of the glutathione peroxidase system 
when examining this relationship, as variation in the GPX1 gene is related to depression risk and significantly 
influences the protective impact of selenium, which is indicative of a gene-environment interaction. 
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Background 

Depression is the most frequently occurring psychiatric 
disorder [1] with 1 in 10 adults reporting depression [2]. 
Depression has been shown to have a wide variety of ad- 
verse health outcomes including poorer physical and social 
functioning, and decreased quality of life [3-5] as well as an 
increase in morbidity and mortality [6-8]. Understanding 
the link between depression and environmental factors that 
may influence its occurrence, symptom presentation, and 
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severity, may have a significant impact on our understan- 
ding of the biological mechanisms underlying depression 
and its treatment. Past research has shown relationships 
between depression and environmental factors, such as 
natural sour gas containing H 2 S [9] and pesticides [10]. 
Another naturally occurring substance that has been 
related to depression is selenium. 

Past research has shown that selenoproteins are essential 
for numerous biological functions [11], many of which are 
thought to be related to depression. Depression is thought 
to be unpinned by multiple pathways, such as inflamma- 
tion, oxidative and nitrosative stress, as well as a decrease 
in antioxidant levels [12-14]. Selenium, in the form of Sele- 
noprotein P, protects the brain against oxidative damage 
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[15] and is involved in the regulation of the inflammatory 
response [11]. Additionally selenium, an antioxidant and an 
essential component in the glutathione peroxidase system, 
protects lipids in polyunsaturated membranes from oxida- 
tive degradation. Depression is often associated with lower 
levels of antioxidants, such as glutathione and GPx, gluta- 
thione peroxidase enzyme [16,17]. These findings suggest 
that selenium may be protective against depression, and 
selenium supplementation may have the potential to alle- 
viate the symptoms of depression. However, research inves- 
tigating the relationship between selenium and depression 
has yielded inconsistent findings. Low selenium concentra- 
tions have been related to increased incidence of many 
human diseases including depression [18-20]. Rayman et al. 
[21] failed to identify any benefit in mood or quality of life 
in a sample of elderly participants. Prenatal selenium sup- 
plementation has been linked to increased serum selenium 
levels as well as decreased risk for postpartum depression 
[22] . Pasco et al. [23] found that lower dietary selenium was 
related to a new diagnosis of major depression. Recently, 
Gao et al. [24] found that higher selenium were associated 
with lower GDS scores; however this effect was no longer 
significant when cognitive functioning was added to the 
model. This finding is interesting because it suggests that 
the link between selenium and depression may be impacted 
by individual factors, such as genetic variants that may 
impact the selenium levels in the body. 

Selenium also plays a key role in proper functioning of 
the glutathione peroxidase system, which has been linked 
to depression. The glutathione peroxidase enzyme (GPx) 
is expressed in lower amounts among depressed patients, 
and GPx activity has been linked to depressed mood and 
autonomic symptoms [25]. Because functional variants 
have been identified within glutathione peroxidase genes, 
genetic variation within the glutathione peroxidase 
pathway is an important consideration when examining 
selenium-health status links. In particular, the Prol98Leu 
polymorphism of the glutathione peroxidase 1 (GPX1) 
gene (rsl050454) has received a great deal of attention re- 
cently in the literature [26] . This polymorphism contains a 
c.599C>T SNP, which is also referred to as rsl050450. 
The single nucleotide polymorphism (SNP) likely affects 
the activity of variant proteins [27]. The T allele of the 
rsl050450 GPX1 SNP has been associated with bladder, 
lung, and breast cancers, intracerebral hemorrhage [28], 
metabolic syndrome [29], and coronary artery disease 
[30]. It is functional in humans and impacts the response 
of GPX1 activity to selenium [31]. Genetic polymorphisms 
of the GPX1 gene have been shown to be related to tumor 
susceptibility [32], blood selenium levels associated 
with supplementation [33], and other health conditions. 
However, to date, we are aware of no studies that have 
examined the impact of GPX1 polymorphisms on the 
selenium- depression link. 



The current study was designed to (a) study the link 
between groundwater selenium levels and depressive 
symptoms, (b) determine if the Prol98Leu GPX1 
directly impact depression risk, and (c) determine if the 
Prol98Leu polymorphisms of the GPX1 gene impact 
these findings. It was hypothesized that higher ground- 
water selenium would be protective against depressive 
symptoms. No specific a priori hypotheses were made 
regarding the Prol98Leu polymorphism given that this 
is the first study to examine the role of this genetic 
variant in depression. 

Methods 

Participants and procedures 

Data were analyzed from 585 participants (183 men and 
402 women) from Project FRONTIER. Project FRONTIER 
is an ongoing study of factors impacting rural health that 
utilizes a community-based participatory research (CBPR) 
approach. The CBPR approach of partnering with commu- 
nity members and organizations is particularly useful when 
working with underserved communities, thus it is very well 
suited to rural health research and has been recommended 
by the National Institute on Environmental Health Sciences 
[34]. All participants signed written informed consents. 
Once participants have completed the consent process, 
s/he undergoes an interview (participant and informant), a 
standardized medical examination, clinical labs, and neuro- 
psychological testing. Inclusion criteria are (a) age 40 and 
above and (b) residing in one of the counties included in 
Project FRONTIER. The counties currently included in 
project FRONTIER are Cochran County, Bailey County, 
and Parmer County, all of which are located on the 
Texas - New Mexico border. Participant recruitment is 
completed by community recruiters/assessors through 
multiple means including brochures/flyers, presenta- 
tions and events, in-person and/or door-to-door solicita- 
tion, as well as snowball recruitment. Our prior work has 
shown that the research sample closely matches the 
demographics of the community at large [34]. Project 
FRONTIER is conducted under an IRB approved proto- 
col. We have published findings utilizing this approach 
with Project FRONTIER elsewhere [35-37]. Project 
FRONTIER is conducted under the approval of Texas 
Tech Health Sciences Center IRB Board, L06-028. 
This archival analysis was conducted from a de-identified 
database under the approval of the University of North 
Texas Health Sciences Center IRB Review Board, protocol 
number 2012-071. 

Measures 

Geriatric depression scale (GDS-30) 

The Geriatric Depression Scale was created specifically 
for screening depressive symptoms among older popula- 
tions [38], and has been widely utilized in both clinical 
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and research settings. The GDS has demonstrated 
adequate internal consistency, test-retest reliability, and 
concurrent validity [39]. A recent factor analytic study 
[40] identified a four- factor solution underlying the 
scale. The Dysphoria factor contains 11 items primarily 
associated with a sad mood; the Meaninglessness factor 
consists of seven items that reflect an appraisal of the 
meaning (or lack thereof) in one's life; the Apathy factor 
is made up of six items that reflect a lack of motivation 
or initiative; finally, the Cognitive Impairment factor 
consists of six items that reflect difficulty and concern 
with cognitive processes. 

Determination of GIS-selenium 

Geographic information system (GIS) analysis is a way of 
displaying and analyzing geographically referenced informa- 
tion that is frequently used to estimate environmental 
exposures [35,41-46]. We used the Environmental Systems 
Research Institute ArcGIS (release 9.2) program to plot 
a point for each of the selenium ground water measure- 
ments from the Texas Water Development Board (TWDB). 
Through inverse distance weighted (IDW) interpolation, 
the ArcGIS software built a three-dimensional surface map 
from a list of points. Each points influence was weighted 
based off its distance to that section of the map, which was 
generated using 12 well measurements from the TWDB 
within the immediate geographic vicinity. Each of the study 
participant's current residential address was geocoded with 
the ArcGIS StreetMap data. Finally, GIS-selenium concen- 
tration was calculated by extracting the estimated selenium 
value from the IDW surface at each resident's location. The 
inter-rater reliability of our GIS approach to calculating our 
groundwater elements is very high (Pearson Correlation co- 
efficient > 0.9). In our prior work, we have found our GIS- 
based estimates to be very similar to concentrations found 
from direct measurement [35]. 

Polymerase chain reactions (PCRs) were utilized to 
genotype the GPX1 (rsl050450) polymorphism. Primers 
were purchased from Integrated DNA Technologies, Inc.: 
FW 5'- CTACGCAGGTACAGC-CGCCGCT-3' and RV 
5'-AAGGTGTTCCTCCCTCG-TAGGT-3'. Amplification 
conditions included: lOOng template DNA, DMSO, water, 
and GoTaq Green Master Mix (3mM MgCl 2 , IX PCR 
buffer, 400 uM each dNTP). The amplicons of the GPX1 
gene were digested with the restriction enzyme Apal with 
0.18 uL per sample, 0.05 uL BSA, and 1 uL lOx buffer 4, 
leaving fragment sizes determined by the genotype of the 
studied individuals (pro/pro = 191 bp; pro/leu = 191, 117, 
and 74 bp; leu/leu = 117 and 74 bp). Digested fragments 
were separated on 1% agarose and visualized by ethidium 
bromide staining. 

In order to examine the link between groundwater sele- 
nium exposures and affective status, linear regression 
models were created with SPSS version 18 using GDS-30 



(total and subfactor) scores as outcome variables and sele- 
nium concentration as the predictor variables. Covariates 
entered into the models were age, gender, education, and 
language of administration. In order to obtain magnitude of 
relations, the amount of variance accounted for by 
selenium concentrations alone is also presented. Logistic 
regression models were created to determine the risk of 
depression, defined as GDS score > 10, that were asso- 
ciated with (a) groundwater selenium concentrations and 
(b) the GPX1 Prol98Leu polymorphism. 

Results 

The mean age and education of the sample was 61.49 (sd = 
12.70; range = 40-96) and 10.76 (sd = 4.37; range = 0-20), 
respectively. Mean GIS-based selenium concentration was 
17.70 ug/L (sd = 10.66; range = 3.96-56.32). Seventy-seven 
percent of the sample was interviewed in English and 33% 
in Spanish; 88% of the sample self-reported their race as 
White, and 38% reported their ethnicity as Mexican 
American. Mean GDS-30 score of the sample was 10.49 
(sd = 5.07; range = 0-26). Demographic characteristics of 
the sample are presented in Table 1. 

Selenium concentrations were found to be significantly 
and negatively related to all GDS scores. Specifically, 
higher selenium levels were significantly related to lower 
scores in the Total GDS-30 scores (|3 = -0.34, p < 0.001), 
explaining 11% of the variance. Higher selenium was 
also related to a decrease in factor scores of Dysphoria 
(|3 = -0.23, p < 0.001), explaining 5% of the variance, 
Meaninglessness (|3 = -0.42, p < 0.001), explaining 17% 
of the variance, and Apathy ((3 = -0.30, p < 0.001), 
explaining 9% of the variance. Higher selenium concen- 
tration was statistically a significant predictor of lower 
scores on the subfactor of Cognitive Impairment (|3 = -0.11, 
p = 0.01), but only accounted for 1% of the variance. The 
results of the regression analysis can be found in Table 2. 
Higher selenium concentration was also associated with a 



Table 1 Demographic characteristics 





Mean (sd) 


Range 


Age 


61.45 (12.70) 


40-96 


Education 


1 0.76 (4.37) 


0-20 


GIS-selenium (ug/L) 


17.70 (10.67) 


3.96-56.32 


GDS Total 


1 0.49 (5.07) 


0-26 


GDS Dysphoria 


3.31 (2.44) 


0-11 


GDS Meaninglessness 


2.24 (1 .43) 


0-7 


GDS Apathy 


2.87 (1 .42) 


0-6 


GDS Cognitive Impairment 


1.87 (1.19) 


0-6 


GPX1 Polymorphism (%) 


CC 46 


CT 48 


TT 


6 
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Table 2 Linear regression results 





B (SE) 


t-score 


p value 


% Variance accounted For 


GDS Total 


-0.16 (0.02) 


-8.21 


0.000 


11% 


GDS Dysphoria 


-0.05 (0.01) 


-5.24 


0.000 


5% 


GDS Meaninglessness 


-0.06 (0.006) 


-10.31 


0.000 


17% 


GDS Apathy 


-0.04 (0.006) 


-7.11 


0.000 


9% 



GDS Cognitive Impairment -0.13(0.005) -2.53 0.012 1% 



significantly lower risk of depression, with B = -0.06 
(SE = 0.01), Wald (df = 1) = 37.25, p < 0.001, and 
OR = 0.94 (95% CI = 0.92-0.96). 

The link between depression and GPX1 polymorp- 
hisms was strongest for CC (homozygous wild type) 
and TT (homozygous mutant type) carriers, with the 
greatest impact being for TT individuals as shown in 
Table 3. The TT polymorphism was the least frequent 
in the sample; however, for TT individuals, selenium 
concentrations accounted for 26% of the variance in 
GDS Total Score, 13% of the variance in Dysphoria, 
37% of the variance in Meaninglessness, and 19% of 
the variance in Apathy. For the CC polymorphism, 
selenium concentrations accounted for 15% of the 
variance in total score, 9% of the variance in Dysphoria, 
22% of the variance in Meaninglessness, and 14% of the 
variance in Apathy. The amount of variance for CT 
polymorphism was relatively small. 

Next, we examined the impact of the GPX1 Prol98Leu 
polymorphism on the findings. First, we sought to 



determine if homozygous carriage of the C allele was asso- 
ciated with depression status (GDS > 10 = depressed), 
compared to T allele carriers (CT plus TT; TT was not 
analyzed separately due to sample size). Logistic regression 
was conducted with depressed status as the outcome vari- 
able and GPX1 polymorphism (CC versus CT) as a 
categorical predictor variable. Sixty-six percent of CC 
homozygotes were depressed compared to 57% for T allele 
carriers; therefore, homozygous carriage of the C allele was 
associated with a significantly increased risk of depression 
(B = 0.50, SE = 0.19, Wald (df = 1) = 6.88, p = 0.009, and 
OR = 1.65 (95% CI = 1.14-2.40). Next, a linear regression 
model was created to examine the gene-dose effect of 
GPX1 C versus T (CC = 2, CT = 1, TT = 0) alleles on 
depression scores (GDS total score). With age, gender, 
education, and test language entered as covariates, the 
increasing presence of the T allele was associated with sig- 
nificantly higher GDS scores (B = 0.94, SE = 0.35, t = 2.69, 
and p = 0.008). Lastly, we re-ran our regression analyses 
examining the link between selenium concentrations and 



Table 3 Linear Regression Results by GPX1 Allele 

B (SE) t score p value Model variance Selenium variance 



CC (n = 221) 












GDS Total 


-0.18 (0.03) 


-6.40 


<0.001 


21% 


15% 


GDS Dysphoria 


-0.05 (0.02) 


-3.24 


0.001 


9% 


4% 


GDS Meaninglessness 


-0.06 (0.007) 


-8.04 


<0.000 


24% 


22% 


GDS Apathy 


-0.05 (0.008) 


-6.10 


<0.000 


18% 


13% 


GDS Cognitive Impairment 


-0.02 (0.007) 


-2.96 


0.003 


7% 


5% 


CT (n = 208) 












GDS Total 


-0.13 (0.04) 


-3.58 


<0.001 


10% 


6% 


GDS Dysphoria 


-0.05 (0.02) 


-2.99 


0.003 


9% 


4% 


GDS Meaninglessness 


-0.05 (0.01) 


-4.68 


<0.001 


12% 


9% 


GDS Apathy 


-0.03 (0.01) 


-2.74 


0.007 


3% 


3% 


GDS Cognitive Impairment 


0.00 (0.008) 


-0.05 


0.96 


2% 


0% 


TT (n = 33) 












GDS Total 


-0.27 (0.07) 


-3.78 


<0.001 


40% 


26% 


GDS Dysphoria 


-0.08 (0.04) 


-2.27 


0.03 


21% 


13% 


GDS Meaninglessness 


-0.11 (0.2) 


-5.11 


<0.000 


53% 


37% 


GDS Apathy 


-0.06 (0.02) 


-2.81 


0.009 


28% 


19% 


GDS Cognitive Impairment 


-0.02 (0.3) 


-0.73 


0.47 


3% 


2% 
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GDS scores (total and subfactors) after stratification 
of the sample on genotype at the GPX1 Prol98Leu 
polymorphism. The link between selenium levels and 
depression was strongest among CC and TT, as compared 
to CT individuals. The TT polymorphism was the least 
frequent in the sample. Of note, in regression analysis of 
groundwater selenium levels on depression risk, models 
that included Prol98Leu genotype along with age, gender, 
education, and language of test administration accounted 
for more variance (70%) than those that did not include 
genotypic status (50%). For CC polymorphism, selenium 
concentrations alone accounted for 15% of the variance in 
total score, 9% of the variance in dysphoria, 22% of the 
variance in meaninglessness, and 14% of the variance in 
apathy, which was more variance accounted for than 
found by covariates. The amount of variance for CT 
polymorphism was relatively small. See Table 3. 

Therefore, these findings point to a gene-environment 
interaction, which is supported by the gene-dose findings. 

Discussions and conclusions 

In the past three decades, selenium has been found to be 
essential to various aspects of health [47]. It has been 
found to play a role in prevention of cancer and cardiovas- 
cular disease as well as in regulating mood [20]. Exposure 
to environmental factors, such as selenium, can also play a 
role in depression and anxiety states. Past research has 
shown that low selenium concentrations were related 
to depression status and other negative affective states 
[18,20]. In addition, studies have indicated that selenium 
may be useful in preventing postpartum depression [22] . 

In the current study, higher groundwater selenium 
levels were associated with lower scores on the Geriatric 
Depression Scale. This finding is consistent with past re- 
search indicating that low selenium status is related to 
increased incidence of depression and other negative 
affective states [18,20]. In addition to the decrease in 
total depression scores, higher selenium levels were spe- 
cifically related to a decrease in symptoms of Dysphoria, 
Meaninglessness, and Apathy, indicating that selenium 
may have an impact on particular types of depressive 
manifestations. It is interesting to note that selenium 
concentrations were only minimally related to the cogni- 
tive impairment subscale, as we have previously shown 
selenium to be selectively protective of memory func- 
tioning [36]. This suggests that, despite the potential role 
of selenium concentrations in preserving cognitive 
function, depressive symptoms related to thoughts of 
cognitive dysfunction appear to be unrelated and may 
be related to the concept of insight. 

This is the first study to examine residential groundwater 
selenium levels and depression. A myriad of factors deter- 
mine the chemical make-up of groundwater in a given area, 



and it is important to be cognizant of the potential health 
impacts from such substances. Research on groundwater 
and its implications on human health can have profound 
implications in the realms of public health and policy. 
Additionally, results obtained by groundwater research 
could prove to be the impetus for the formation or refor- 
mation of population-based prevention strategies concer- 
ning groundwater safety and exposure. This study was 
conducted in a rural agricultural area where higher concen- 
trations of groundwater selenium are typically found. 
Aligning with this, the mean selenium level in this study 
was 17.70 ug/L, whereas the national average is 10 ug/L 
[48]. 

The protective impact of selenium on mental health fits 
with its biological activity and is supported by both animal 
model and human research. Selenium, a main constituent 
of neuronal selenoproteins, may be particularly important 
for sustaining healthy cerebral function [47] through its key 
role in protection against oxidative damage [15,20,49]. Oxi- 
dative damage is a key biological mechanism in depression 
pathophysiology, and many anti-depressants, like fluoxetine 
[50] and venlafaxine [51], reduce oxidative damage. An- 
other possible mechanism of action is through selenium's 
impact on brain-derived neurotrophic factor (BDNF), which 
has been linked to depression and found to be decreased in 
selenium deficient rats [52]. Additionally, depression is also 
associated with lower antioxidant levels. Research has sug- 
gested that antioxidant deficiency found in depression may 
predispose individuals to greater inflammatory and oxida- 
tive and nitrosative responses [13], which may contribute to 
the protective effects of selenium against depression. 

Ours is the first study to directly examine the link 
between the GPX1 Prol98Leu polymorphism and depres- 
sive risk as well as its impact on the selenium-depression 
link. Our findings point to gene-environment interaction, 
which is supported by the gene-dose impact. Specifically, 
the C allele of the GPX1 polymorphism was associated 
with significantly increased risk of depression in a dose- 
dependent manner (i.e. CC > CT > TT). Additionally, the 
impact of environmental selenium concentrations on de- 
pression scores was stronger for the CC than CT group. 
The TT group showed a substantial impact in the expected 
direction, but due to the small sample number of TT 
homozygotes, these results must be replicated. It is also im- 
portant to note the differences in effect sizes for selenium 
on depression observed in each genotypic group. Among 
CC homozygotes, the majority of variance accounted for in 
GDS scores was from selenium concentrations with age, 
gender, education, and language of test administration con- 
tributing less to GDS scores. Among the TT group, the 
amount of variance accounted for by selenium concentra- 
tions alone was substantial; however, covariates in the 
model also accounted for a larger portion of the variance. 
Future work is needed to delve more thoroughly into this 
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gene-environment interaction as well as the dose-response 
impact of the C allele of the GPX1 polymorphism on risk 
for depression. Additional work is also required to examine 
other potential genes of interest. For example, Galecki and 
colleagues [53] identified polymorphisms of the NOS2A 
gene, which are important in the nitric oxide pathways, as 
related to recurrent depression. Given that selenium (i.e. 
selenoproteins) is involved in the nitric oxide pathway as 
well, NOS2A and other genes relevant to pathways linked 
to selenium should be examined. It is likely that more 
sophisticated bioinformatic analyses will be necessary to 
disentangle these multi-level gene-(biological mechanism) 
- environmental exposure - depression interactions. 

These results point towards possible pharmacogenetic 
implications of the Prol98Leu GPX1 polymorphism. Spe- 
cifically, in our sample, CC homozygotes had significantly 
increased risk for depression; however, selenium was also 
most powerfully protective against depression and depres- 
sive symptoms among those particular carriers. Therefore, 
future research looking at the therapeutic impact of se- 
lenium on depression should examine this genotype. 
Additionally, previously conducted clinical trials of anti- 
depressant medications with stored DNA (or previously 
conducted genetic analyses) can determine if Prol98Leu 
GPX1 genotype predicts optimal treatment response. 
Lastly, the Prol98Leu GPX1 polymorphism may be useful 
for identifying those individuals most likely to benefit 
from anti-depressant medications in new clinical trials. 

There are limitations to the current study. One is the 
GIS-based estimate of selenium concentration rather than 
direct measurement. However, our prior work, as well as 
that of others, has shown that GIS-based methods are valid 
ways for estimating environmental exposures. The current 
findings should be replicated using direct selenium meas- 
urement. The cross-sectional nature of the study is another 
limitation; however, our prior work has shown selenium to 
be preventative of cognitive decline prospectively in a sub- 
set of this cohort [30]. Future studies will examine the 
protective impact of selenium concentrations on both 
depression and cognitive decline as future waves become 
available. The current study also has several strengths. 
Specifically, this is the first study to examine this link 
in such detail (i.e. depressive symptoms rather than 
only global scores) as well as look at this topic among 
rural- dwelling adults and elders. This is also the first 
study to examine the impact of genetic variation within 
GPX1 on depression risk as well as the selenium- 
depression link. Taken together, these results point to the 
need for continued work in the area, as selenium may 
provide a therapeutic target for preventative strategies 
aimed at reducing depressive incidence rates in vulnerable 
populations. Additional work is also needed to determine 
the pharmacogenetic implications of the Prol98Leu GPX1 
polymorphism in depression therapies. 
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